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Abstract 

For a high temperature non-Abelian plasma, we reformulate the hard thermal 
loop approximation as an effective classical thermal field theory for the soft modes. 
The effective theory is written in local Hamiltonian form, and the thermal partition 
function is explicitly constructed. It involves an ultraviolet cutoff which separates 
between hard and soft degrees of freedom in a gauge-invariant way, together with 
counterterms which cancel the cutoff dependence in the soft correlation functions. 
The effective theory is well suited for numerical studies of the non-perturbative 
dynamics in real time, in particular, for the computation of the baryon number 
violation rate at high temperature. 
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The violation of the baryon number in the high-temperature, symmetric phase of 
the electroweak theory is an important example of a physical process which is sensitive to 
the non-perturbative real-time dynamics of hot gauge theories [[I]. Unlike static charac- 
teristics, like the free energy, which can be computed on the lattice, the non-perturbative 
evolution in real-time cannot be studied through the standard lattice simulations for- 
mulated in imaginary-time. However, it has long been recognized that a fully quantum 
calculation is actually not necessary Q: the non-perturbative phenomena are associated 
with long wavelength[] (A ^ l/g 2 T) magnetic fields, which, because of Bose enhancement, 
have large occupation numbers, 

N °W = jwz^-is for E<<T ' (1) 

and should therefore exhibit a classical behaviour. Based on this observation, there has 
been attempts to compute the baryon number violation rate T through lattice simulations 
of the classical thermal Yang-Mills theory |2|, E|. 

However, it has been recently observed j|, |5| that the baryon violating processes are 
actually sensitive to the hard thermal modes with momenta ~ T, for which the classical 
approximation is well-known to fail: the hard modes cause the damping of the soft field 
configurations, an effect which is predicted to reduce T by a factor of g 2 as compared 
with its classical estimate in Refs. [[[], [2|] . In order to verify this prediction and eventually 
compute T, one has to properly take into account the effects of the hard modes on the 
dynamics of the soft fields. To leading order in g, these effects are encompassed by the 
so-called "hard thermal loops" (HTL) [6-11], which are non-local one-loop corrections to 
the soft (k ^ gT <C T) field propagator and vertices due to the hard (k ~ T) thermal 
modes. 



It has been first suggested in Ref. |12[ to use the local version of the HTL effective 
theory [8-11] in classical lattice simulations in order to compute the baryon number 
violation rate. However, in order to transpose this idea into practice, one first needs a 
precise, and gauge-invariant, separation between hard and soft degrees of freedom (to 
avoid overcounting, and to provide an ultraviolet cutoff to the effective theory for the soft 
modes). Loosely speaking, this requires an intermediate scale fi, with gT C/iCT, which 
should act as an infrared (IR) cutoff for the hard modes and as an ultraviolet (UV) cutoff 

denotes the temperature, assumed to be large enough for the coupling constant g(T) to be small: 

9 < 1. 
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for the soft ones, and which should cancel in the calculation of physical quantities. But the 
practical implementation of such a separation of scales meets with technical difficulties: 
(i) In the hard sector, one cannot simply introduce /i as an infrared cutoff in the one-loop 
diagrams for the HTL's since this would break gauge symmetry |Tj|. (ii) In the soft sector, 
one cannot use the lattice spacing to provide the necessary upper cutoff ~ fi: indeed, a 
finite (and relatively large: a ~ ^> T _1 ) lattice spacing introduces lattice artifacts 
which make impossible the matching with the hard sector 0, 

Another problem is the proper definition of classical thermal expectation values 
within the effective theory. In order to compute such expectation values, one has to 
solve first the equations of motion for given initial conditions, and then average over the 
classical phase space with the Boltzmann weight exp(—(3H). Since the local formulations 
of the HTL theory which are currently available || |9|, |ll|] are not in canonical form, it is a 
non-trivial task to identify the independent degrees of freedom and construct the classical 
phase space. 

The purpose of this Letter is to give an explicit solution to the above mentioned 
problems by constructing a local effective theory for the soft modes with a /i-dependent 
Hamiltonian. The classical partition function will involve an explicit UV cutoff, chosen 
so as to cancel the //-dependence of the Hamiltonian in the calculation of IR-sensitive 
correlation functions, like the baryon number violation rate T. 

Our construction relies in an essential way on the local formulation of the HTL 
theory presented in Ref. ||, which we briefly review now. It involves a set of coupled 
equations for the soft fields and their induced current, namely eqs. (@)-(fD below. The 
soft fields A*(x) satisfy the Yang-Mills equations with an induced current in the right 
hand side: 

(D u F"») a = ft, (2) 

where Z> M = <9 M + ig[A^ ■ ], A^ = A a ^T a , and = [D^, D u ]/(ig). (The generators of 
the colour group in the adjoint representation are denoted by T a ; they satisfy [T a ,T b ] = 
i: fabc T c and Tr ( T « T b) = C A 5 ab , with C A = N for SU(iV).) The induced current j M = j*T a 
is related to the colour fluctuations of the hard thermal modes: 

f a {x) = 2gC A J ^L v »5N a (k,x). (3) 

In this equation, 5N(k, x) = 5N a (k, x) T a is a phase-space colour density matrix for hard 
gluons (|k| ~ T), which describes long- wavelength colour correlations as induced by the 
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soft fields A". Furthermore, = (1, v) and v = k/k is the velocity of the hard particle 
(k = |k|, and |v| = 1). The system is closed by the kinetic equation for the density 
matrix, which is a non-Abelian generalization of the Vlasov equation: 

(v-D x )5N(k,x) = -gy.E(x)^-, (4) 

where E l a = and N Q (k) = l/(e l3h — 1). The dynamics described by the above equations 
is gauge invariant, and the current j" is covariantly conserved: D^j^ = 0. 

From eq. (^j), we note that the v and fc-dependence can be factorized in <5iV a (k, x) 
by writing: 

6N a (k, x) = -gW a (x, v) (dN /dk). (5) 
The new functions W a (x, v) satisfy the equation: 

(v-D x )W(x,v) = v-E(4 (6) 

which is independent of k since the hard particles move at the speed of light: |v| = 1. 

Eqs. above provide a local description of the soft field dynamics in the HTL 

approximation. In order to use these equations for classical thermal calculations, one needs 
to (i) introduce an IR cutoff [i in the hard sector, (ii) perform a Hamiltonian analysis (to 
identify the independent degrees of freedom and the corresponding Hamiltonian), (iii) 
write down the classical partition function, and (iv) supply the effective theory with an 
UV cutoff ~ /i. We shall address these problems in this order: 

(i) By inspection of eqs. (@)-@, it is quite obvious how to introduce the intermediate 
scale fi: since k is the momentum carried by the hard particles, it is sufficient to integrate 
in eq. (|3D with a lower cutoff equal to fi. With eq. (Q), the radial integration in eq. (0) 
can be worked out, with the result 

fa(x) = m 2 Mj^v»W a (x,v), (7) 

where the angular integral / dfl runs over the unit sphere spanned by v, and 

The quantity = = 0) is the physical Debye mass to leading order in g || 0, 

For gT <C \i <T, mn{p) is the hard sector contribution to mo- 
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In contrast to the usual one-loop calculations [E2|, the above implementation of /i 
has preserved gauge symmetry automatically: indeed, the kinetic equation ([|) is gauge 
covariant for any value of k, so that the /z-dependent current in eq. (|7|) is covariantly 
conserved. 

(ii) A Hamiltonian analysis of the HTL theory has been given by Nair |J, in terms 
of some new auxiliary fields. Here, we shall rather follow Refs. [11] and propose a simpler 



Hamiltonian formulation which involves the fields W a (x, v) introduced above. In the gauge 
Aq = 0, the independent degrees of freedom are E£, A\ and W a , and the corresponding 
equations of motion follow from eqs. @, (H) and (0) above: 



E? = -d Al 



-d Et + e ijk (D j B k ) a = m 2 Mj — Vi W a (x,v), 

(d + v-B) ab W b = v-E a , (9) 

together with Gauss' law which in this gauge must be imposed as a constraint: 

(D-E) a + m^)/|V(*,v) = 0. (10) 

Note that eqs. (|9|) are not in canonical form: this is already obvious from the fact that 
we have an odd number of equations. Accordingly, it is not a priori clear that these 
equations are Hamiltonian in any sense|. It has been shown in Ref. JXTJ] that eqs. (§) are 



conservative: the associated, conserved energy can be computed in any gauge as: 

H = ^|d 3 x|E a -E a + B a B a + m 2 M J ^ W a (x, v) W a (x, v)| . (11) 

Remarkably, we show now that, in the gauge Aq = 0, the functional (|ll|) also acts as a 
Hamiltonian, that generator of the time evolution. To this aim, we introduce the 

following Poisson brackets^ (see also Ref. 0) : 

{£f(x),4(y)} = -5<%<5(3)( x -y), 
{^(x),H/ 6 (y,v)} = v t 5 ab 5^^-y), 
m ^{r(x,v),#(y,v')} = (gr bc W c + (vB x ) ab )5®(x-y)5(vy). (12) 



* Recall that the Hamiltonian structure is a non-trivial issue already for the Abelian Maxwell- Vlasov 
equations jL4j, of which eqs. (^|) can be seen as a non- Abelian generalization. 

§More precisely, these are generalized Lie-Poisson brackets, according to the terminology in Ref. 
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Here, <5(v, v') is the delta function on the unit sphere, normalized such that 

|^5(v,v')/(v) = /(v'), (13) 

and all the other Poisson brackets are assumed to vanish. We also assume standard 
properties for such brackets, namely antisymmetry, bilinearity and Leibniz identity. It is 
then straightforward to verify that (a) the Poisson brackets (|l^) satisfy the Jacobi identity 
(as necessary for consistency) and (b) the equations of motion (|9[) follow as canonical 
equations for the Hamiltonian (0). For instance, d^W 0, = {H, W a }, and similarly for E" 1 
and A". 

Note that the effective theory in eqs. ©-(0) involves the infrared cutoff (and 
also the temperature T) only through a single mass parameter, namely the "hard" Debye 
mass m 2 H {ji) of eq. (|]). 

(iii) We are now in position to construct (generally time- dependent) thermal ex- 
pectation values within the classical field theory defined by eqs. (^J)-(|ll|). The thermal 
phase-space is defined by the initial conditions for eqs. (J^j, and the canonical weight is 
given by the effective Hamiltonian (|11|). Thus, the thermal correlation functions of the 
fields A % a can be obtained from the following generating functional: 

Z d [J?] = Jv£?VA?VW a 5(g a )expi [ -f3H + |d 4 xJf(x)4 a (x)| , (14) 

where A\(x) is the solution to eqs. (H) with the initial conditions {S^A®, W a } (that is, 
E^(to, x) = (x), etc., with arbitrary to), and H is expressed in terms of the initial fields. 
Since the dynamics is gauge- invariant, it is sufficient to enforce Gauss' law at t = to: 

g a = (v^r + m 2 M J^-W a = 0. (15) 

The only subtle point in eq. fll"4]) is the definition of the measure in the phase-space[]: 
since we are not using canonical variables, we still have to verify that the naive measure 
T>£^T>AiT>yV a is indeed the correct one. A necessary condition is that this measure be 
invariant under the time evolution described by eqs. (||), so that Z c i[J] be independent of 
to, as it should. This condition can be most easily verified by considering an infinitesimal 
time evolution of the form $ Q — ■> <$>' a = $ Q + {if, $ a }dt, where $ Q refers to any of the field 
variables {£", Af, W a }. Then, by using eqs. (0) — or, equivalently, the Poisson brackets 
am grateful to Tanmoy Bhattacharya for an illuminating discussion on this point. 



5 



(72) — , it is straightforward to verify that the Jacobian for this transformation is equal 
to one, to linear order in dt : 



J 



5(S u Ai,W) 



i + o((dtn (i6) 



As a further check, one can verify that eq. ([14]) reduces to standard results in some 
particular cases. For instance, for J? = this equation yields the result expected from 



dimensional reduction [15|, that is, 

Z d = J VA%VAl exp |-| J d 3 x (B?B? + (ViAoWViAo)* + m 2 MA a A%) j , (17) 

where the Aq components of the gauge fields have been reintroduced as Lagrange mul- 
tipliers to enforce Gauss' law, and the functional integrals over Sf and W a have been 
explicitly performed. Converserly, the general formula fll4] ) can be seen as a general- 
ization of the dimensional reduction method to include dynamical (i.e., time-dependent) 
phenomena. Another check is provided by the Abelian limit, where eq. ( ]T4]) yields, after 
a straightforward calculation : 

Z d [Ji] = expj-i j ' d 4 x j ' d A y J i (x)*D ij (x - y)Jj(y) 

*D l3 {x-y) = J-^e-^-y^pMNMo), (18) 

where *Pij(q) is the magnetic photon spectral density in the HTL approximation || and 
N c i(q ) = T/q is the classical thermal distribution function, which coincides with the low 
energy limit of the quantum distribution (cf. eq. (|1|)). In the second line of eq. ( p~8"D we 
recognize, as expected, the classical limit of the soft photon 2-point function in the HTL 
approximation. 

(iv) The last step toward a well-defined classical effective theory is to supply the 
partition function ( |14) with an ultraviolet cutoff A c / ~ fi, chosen so as to cancel — in the 
calculation of the soft correlation functions — the explicit /i-dependence of the effective 
Hamiltonian (|IID (a cancellation to be subsequently referred to as matching). Since 
the effective theory is ultimately intended for non-perturbative calculations, it will be 
convenient to chose a regularization method which can be also implemented on a lattice. 
This cannot be the lattice spacing itself: indeed, if we choose to work with a finite lattice 
spacing a ~ 1///, then we break rotational and dilatation symmetry, and the UV structure 
of the lattice theory gets so complicated that the matching cannot be performed anymore 
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T2|, |l"3f. Lattice artifacts can be eliminated only by taking the continuum limit a — > 



which requires a to be independent of /i. 

The strategy that we propose here is thus the following: the effective theory will 
be formulated as a cutoff theory in the continuum, and the matching will be performed 
in the continuum, rotationally-invariant theory. Then, for computational purposes, the 
resulting cutoff theory must be put on a lattice with small lattice spacing a <C l/A^ |. 
Because of the explicit UV cutoff A c /, the continuum limit a — > is well-defined. 

Following Ref. |13| |, we introduce a smooth UV cutoff in the continuum theory by 
replacing, in the effective Hamiltonian ([TT|) , 

D 2 

A 1 



TiB i B i — ► Tr£7 ( J B\ (19) 



.2 



where f(z) = 1+z and the trace refers to color indices. Besides being gauge- invariant, the 
regularization prescription in eq. (|19|) has also the advantage that it can be carried out on 



the lattice (by using improved lattice Hamiltonians |13|, [U||). However, this prescription 
breaks down dilatation symmetry and, as a consequence, the matching cannot be per- 
formed for all the soft correlation functions (see below), but only for the non-perturbative 
quantities which are infrared sensitive, like the baryon number violation rate T. Let us 
explain this in more detail: 

For matching purposes, we need the A c /-dependent corrections to the soft correlation 
functions to one-loop order in the effective theory. Rather than computing loop diagrams, 
it is more convenient to rely on kinetic theory to describe the interactions between the 
genuinely soft fields, with momenta k ^ gT, and the relatively "hard" classical modes, 
with momenta k ^ A c /. The relevant kinetic equations can be derived in the same way 
||, and look similarly, to our previous equations ([|H(|])- The only differences refer to the 
replacement of N (k) by N c i(Ek) = T/Ek, and of the unit vector v by the group velocity 
Vfc = Vfc-Efc- (These prescriptions can be readily verified by inspecting the derivation 
of the kinetic equations in Ref. H; they have been also justified by a diagrammatic 
analysis in Ref. |13[.) Here, E\ = k 2 f (k 2 / K 2 cl ) is the dispersion equation for the relatively 
"hard" (k ^ A c z) classical excitations, for which the HTL corrections are relatively small 
(since win ~ gT A c z) and can be neglected to the order of interest. It then follows 



"Note that a similar strategy has been recently proposed by Arnold, in the context of purely classical 
Yang-Mills theory jL3|; of course, the matching was not an issue in Ref. [i3|] , which was rather concerned 
with improving the rotational symmetry in classical lattice simulations. 
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that the colour current due to the "hard" classical modes — and which summarizes the 
A c ;-dependence of the effective theory to one loop order — has the form (compare to 
eqs. 0-0): 

#•<*) = -2 9 ^/|^ Wl ,v t ), (20) 

with the functions W a (x, v k ) satisfying (in the temporal gauge A = 0): 

(d + v k -B) ab W b (x,v k ) = v fc -E a (x). (21) 

An important difference with respect to eq. (0) is that the classical "hard" excitations do 
not move at the speed of light, but rather with a fc-dependent velocity v&. Accordingly, 
the radial and angular integrations in eq. ( ^0|) cannot be disentangled anymore, and the 
current jg a will not be characterized, in general, by a single mass scale (in contrast to the 
HTL current in eq. (|7D), which prevents us from performing a full matching. 

However, as we show now, the matching can still be done in the calculation of 
non-perturbative quantities which are infrared sensitive. It is indeed well-known (see, 



e.g., Refs. 0, ^, [13], [T7|]) that the field configurations which are responsible for the non- 
perturbative phenomena (and also for the IR divergences of the perturbation theory) are 
very soft (k = |k| ~ g 2 T) magnetic fields of almost zero frequency: k ^ g 4 T <C k. 
Indeed, these are the only configurations which are not screened at the scale gT by the 
HTL's. For such fields, time derivatives are suppressed with respect to spatial gradients, 



and eq. (pip reduces to: 

(v-B)W(x,v k ) = vE(x), (22) 

where we have been able to simplify one factor of |v&|, so that v is an unit vector, as in 
eq. ([|). Eqs. (|22|) shows that on the relevant, non-perturbative field configurations, the 
function W(x, v k ) = W(x, v) is independent of the radial momentum k. Then, the radial 
integral in eq. (|20|) factorizes and yields (for the relevant, magnetic piece of the current): 

f s a (x) ~ m|(A d ) J^v* W a (x, v), (23) 

with (compare to eq. (|)): 



m s (A d ) = — / dkk v fc 

n z Jo 

2g 2 C A 



d 



IT 2 JO 



dE k 

POD 

/ dkkN d (E k ) = Kg 2 C A TK cl . (24) 
Jo 



k 

oo 



The precise value of the numerical coefficient k can be found in eq. (4.8) of Ref. [JT3[ . 

Eq. ( |23|) shows that, to one- loop order in the effective theory, the A d - dependent 
corrections to the amplitudes involving quasistatic and soft (k <C k ~ g 2 T) magnetic 
fields are characterized by a single mass scale, namely m|(A c ;). This is similar to the HTL 
current in eq. (|7|), so it is now possible to perform the matching by requiring m|(A c ;), 
eq. (fj4|), to cancel the /i-dependent piece of m 2 H (fi), eq. (H). This is achieved by choosing 
[i = n 2 K,A c i : with this matching condition, the IR-sensitive quantities computed in the 
effective theory (|^)-(|T4]) with the UV regularization (|i~9"D come out independent of \x and 



A d , for fi in a large range of values: gT <C \i < T. 

In particular, the effective theory thus defined provides a /i-independent value for 
the baryon number violation rate T, to be ultimately computed on the lattice. The lattice 
implementation of the effective theory (which requires a lattice version of the new fields 
W a (x,v) and of the corresponding equations of motion) is by itself a non-trivial issue, 
which remains beyond the scope of the present work. 

Let us finally note a different proposal |18[] for including the HTL's, which is to treat 



the hard degrees of freedom as classical coloured particles [ID]. This has been recently 



implemented in lattice simulations with results which seem to confirm the predictions 
in Ref. || . By comparaison, the method that we have proposed here, besides being derived 
from first principles, has also the advantages to give the hard modes the correct quantum 
statistics, to involve no free parameter, and to perform a precise matching between hard 
and soft degrees of freedom, thus allowing for the continuum limit to be taken in lattice 
simulations. 

To conclude, we have provided an effective classical thermal field theory for the soft 
modes which includes the hard modes in the HTL approximation and which is well-suited 
for numerical studies of the real-time non-perturbative dynamics. Important applications 
include the high-T anomalous baryon number violation, the dynamics of the electroweak 
phase transition (which requires adding the Higgs field to the above theory), and non- 
perturbative properties of a hot quark-gluon plasma. 

Acknowledgements: During the elaboration of this paper, I have benefited from discus- 
sions and useful remarks from a number of people. It is a pleasure to thank J. Ambj0rn, 
T. Bhattacharya, J. P. Blaizot, S. Habib, A. Krasnitz, L. McLerran, E. Mottola, J.Y. 
Ollitrault, M.E. Shaposhnikov, A. Smilga and N. Turok. 



9 



References 



[2 

[3 

[4] 
[5] 
[6 

[7 



[9 
[10 

[11 
[12 

[13 

[14 



For a recent review, see V.A. Rubakov and M.E. Shaposhnikov, |hep-ph/9603208 
Usp. Fiz. Nauk. 166 (1996) 493. 

J. Ambj0rn and A. Krasnitz, Phys. Lett. B362 (1995) 97; ibid. Nucl. Phys. B506 
(1997) 387. 

G. D. Moore and N. G. Turok, Phys. Rev. D 55 (1997) 6538; ibid. Phys. Rev. D56 
(1997) 6533. 

P. Arnold, D. Son, and L. Yaffe, Phys. Rev. D55 (1997) 6264. 
P. Huet and D. Son, Phys. Lett. B393 (1997) 94; D. Son, |hep-ph/9707351 



R.D. Pisarski, Phys. Rev. Lett. 63 (1989) 1129; E. Braaten and R.D. Pisarski, Nucl. 
Phys. B337 (1990) 569. 

J. Frenkel and J.C. Taylor, Nucl. Phys. B334 (1990) 199; J.C. Taylor and S.M.H. 
Wong, ibid. B346 (1990) 115; R. Efraty and V.P Nair, Phys. Rev. Lett. 68 (1992) 
2891; R. Jackiw and V.P. Nair, Phys. Rev. D48 (1993) 4991. 

J. P. Blaizot and E. Iancu, Phys. Rev. Lett. 70 (1993) 3376; Nucl. Phys. B417 (1994) 
608. 

V.P. Nair, Phys. Rev. D48 (1993) 3432; ibid. D50 (1994) 4201. 

PF. Kelly, Q. Liu, C. Lucchesi and C. Manuel, Phys. Rev. Lett. 72 (1994) 3461. 

J. P. Blaizot and E. Iancu, Nucl. Phys. B421 (1994) 565; ibid. B434 (1995) 662. 

D. Bodeker, L. McLerran, and A. Smilga, Phys. Rev. D52 (1995) 4675. 

P. Arnold, Phys. Rev. D55 (1997) 7781. 

J.E. Marsden and T.S. Ratiu, Introduction to Mechanics and Symmetry: A Basic 
Exposition of Classical Mechanical Systems (Springer Verlag, 1995). 



[15] K. Farakos, K. Kajantie, K. Rummukainen and M.E. Shaposhnikov, Nucl. Phys. 
B425 (1994) 67; K. Kajantie, M. Laine, K. Rummukainen and M.E. Shaposhnikov, 
Nucl. Phys. B458 (1996) 90; E. Braaten, Phys. Rev. Lett. 74 (1995) 2164. 

10 



[16] G.D. Moore, Nucl. Phys. B480 (1996) 689. 

[17] J. P. Blaizot and E. Iancu, Phys. Rev. Lett. 76 (1996) 3080; Phys. Rev. D56 (1997) 
7877. 

[18] C. Hu and B. Muller, Phys. Lett. B409 (1997) 377. 
[19] G.D. Moore, C. Hu and B. Muller, |hep-ph/9710436| . 



11 



